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We propose a simple yet efficient method for generating in-plane hollow beams with a nearly-full circular 
light shell without the contribution of backward propagating waves. The method relies on modulating the 
phase in the near field of a centro-symmetric optical wavefront, such as that from a high-numerical-aperture 
focused wave field. We illustrate how beam acceleration may be carried out by using an ultranarrow non-flat 
meta-surface formed by engineered plasmonic nanoslits. A mirror-symmetric, with respect to the optical axis, 
circular caustic surface is numerically demonstrated that can be used as an optical bottle. 
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I. INTRODUCTION 

Optical bottle beams that contain an intensity null 
point generated considerable interest due to their po¬ 
tential application in optical confinement of laser-cooled 
atoms and in optical tweezers systems^ - — . This bottling 
effect appears when plane waves in the Fourier expansion 
of beams interfere destructively to generate a bottle or 
a null intensity point at the center of the beam. Con¬ 
ventional hollow beams used in optical tweezers and fol¬ 
lowing the abovementioned basic concept are Laguerre- 
Gaussian beams. More recently, an optical scheme for 
realizing a supersized light capsule was presented theo¬ 
retically and experimentally in air, where the volume of 
darkness inside the light capsule can be changed in a wide 
ranged. In this case, a binary phase mask was used to 
create such a tunable antiresolution region. 

Alternatively, wave fields exhibiting a peak intensity of 
curvilinear trajectory that additionally constitutes the 
closed boundary of a geometrical shadow, instead of a 
null point, have been proposed to be used as optical 
bottles^. A well-known example of this kind of wave 
fields are the so-called Airy beams. Unlike ordinary opti¬ 
cal wavefronts, Airy beams transversely accelerate, pre¬ 
serving their intensity profile throughout propagation but 
the profiles features follow a curved path^. By appro¬ 
priately superimposing mirror-symmetric Airy beams, 
the field amplitude oscillates outward of a dark closed 
region^. These theoretical predictions were subsequently 
verified by experimental observations^. The resulting 
caustic surface that abruptly emerges and subsequently 
autofocuses was used as an optical bottle to trap and 
manipulate dielectric microparticlesii - — . 
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The intensity peak of an Airy beam follows a parabolic 
trajectory much like the ballistics of projectiles. The 
shape of the optical bottle produced by this sort of parax¬ 
ial accelerating beam is also governed by the parabolic 
caustic curve. This fundamental limit may be overcome 
since Fourier-space generation of arbitrary accelerating 
beams allow the propagation along controlled curvilin¬ 
ear trajectories^ - — . For instance, nonparaxial accel¬ 
erating waves in two dimensions, whose intensity pro¬ 
files are roughly preserved within a range of propaga¬ 
tion distances, and whose maxima follow semicircular 
paths, were recently propose d 17 ! 18 . Such a beam with 
Bessel signature should be able to bend from a launch 
angle of +90° —perpendicular to the original direction 
of propagation— all the way to an angle of zero, and 
immediately accelerating to propagate at angles close to 
—9Q o19 ! 20 . The wave field will propagate backward in the 
remainder half circular trajectory, a fact that in practice 
limits the formation of light capsules of circular symme¬ 
try. 


In this study we extend these ideas to synthesize an op¬ 
tical bottle exhibiting all-inclusive circular symmetry, in¬ 
tegrating forward-only propagating waves. To construct 
these waves we combine two concentric incomplete Bessel 
beams by applying the concept of symmetrization in a no¬ 
tably dissimilar way of that discussed in Ref^i. For that 
purpose we follow a method previously shown in Ref/ 22 
that potentially is appropriate for using in a compact¬ 
sized system, which makes use of a non-planar meta¬ 
surface sustained by dispersion localities. Our multilay¬ 
ered metal-dielectric nanostructure is piecewise periodic 
in the angular coordinate, thus enabling to transform a 
high-aperture focused beam into a mirror-symmetric ac¬ 
celerating beam nearly enclosing a circular region. 
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II. THEORETICAL PRELIMINARIES 

Our study is based on the efficient transformation 
of a converging wave field into a mirror-symmetric ac¬ 
celerating beam with a prescribed circular trajectory. 
Such beam shaping is achieved by a meta-surface, which 
diffracts the impinging circular wavefront to be modu¬ 
lated appropriately. For that purpose, let us consider 
the diffraction of a monochromatic converging wave by 
a diffracting optical element. In this paper, for the 
sake of clarity, we will consider two-dimensional waves 
propagating in the xy plane. However, a generaliza¬ 
tion in three dimensions may be carried out straightfor¬ 
wardly. The polarization of the incident focused field 
and the subsequently diffracted field will be transverse 
magnetic (TM), in such a way that the magnetic field 
with time-domain frequency uj can be set as H (x,y,t) = 
H(x,y)exp(—iu)t)z, where z is the unitary vector point¬ 
ing along the z axis. In addition, this will be of help in 
order to excite surface plasmon resonances in the elemen¬ 
tary nano-slits taking part of the designed meta-surface. 



FIG. 1. Schematic arrangement representing the diffracted 
converging wave of focal point F as evaluated from the Debye 
diffraction Eq. ©• The emerging wave field propagates in the 
xy plane, and the meta-surface induces a wavefront deviation 
from its original shape (reference cylinder) by w{0). A beam 
carrying a given acceleration leads to wave localization along 
an incomplete circumference with center at focus F. 

According to the Debye diffraction theory, the wave 
field in the focal region of the shaped converging beam 
can be obtained by means of the following diffraction 
integra l 23 ! 24 



— 7r/2 


where r = (x, y), k = uo/c is the wave number, and H s 
is the magnetic field over the reference cylinder. As it is 
shown in Fig. [TJ R is the radius of the diffracted cylin¬ 
drical wave as taken over the reference cylinder. Finally, 
q = (sin 0, cos 0) is a unit vector pointing from the fo¬ 
cus F in the direction of a given point P on the wave- 
front. The scattered wave field H S (Q) will be expressed 
by means of a real and positive term, A(0), and a phase 


only term exp [iw(0)\, giving 

H s (6) = A{0) exp [iw(0)\, (2) 

The apodization function A{6) takes into account the 
truncation and local attenuation of the converging field 
after passing through the diffracting optical element, and 
w(6) denotes any deviation of the diffracted wavefront 
from the reference cylinder, which commonly is inter¬ 
preted as monochromatic aberration. 

The apodization function in our simulations is a super- 
Gaussian function instead of a customary flat-top profile, 
which can be written as, 

A(0) = Ao exp[—(#/f2) 6 ], (3) 

of semi-aperture angle £2, enabling to minimize edge ef¬ 
fects. We should mention that the limits of integra¬ 
tion in Eq. 0 including forward-only propagating waves 
makes the diffracted fields practically unaltered for semi¬ 
apertures 12 = 7r/2 and higher values. Furthermore we 
will consider a linear phase modulation, which induces 
the acceleration of the diffracted field around the focal 
point F 22 . By neglecting a constant term of Taylor ex¬ 
pansion of w(Q), which has no significant contribution in 
the diffraction integral Eq. (pQ) , and also neglecting higher 
orders of this expansion, we can rewrite the phase term 
in Eq. © as exp(im6). In fact, the parameter m governs 
the curvilinear trajectory of the focused field around the 
geometrical focus F, as demonstrated below. Particu¬ 
larly, the distance between the focal point F and where 
the beam is localized, is determined by m. It can be 
used for controlling the acceleration of a focused beam. 
Before the diffraction broadening overcomes, the beam 
propagates along a circular trajectory whose radius can 
be estimated as r rn = |m|/fe2£. 

Analytical solution of the Debye integral Eq. © can 
be achieved provided that 12 oo and including back¬ 
ward propagating waves. By considering the parameter 
m as an integer number, the magnetic field H( r) finally 
yields an expression that is proportional to the Bessel 
function, namely J m (fcr), where r = |r| 22 . Disregard¬ 
ing backward propagation, the resulting wave resembles 
a sector Bessel field centered on the positive (negative) 
x axis for m > 0 {m < 0). Because of that, this type 
of accelerating beams can be called as incomplete Bessel 
beamsi^. In Fig. [2ji and (b) we show the intensity dis¬ 
tribution | if | 2 of a Bessel field of order m = 10 and 
the associated incomplete Bessel beam of semi-aperture 
12 = 7r/2. By changing the sign of the parameter m, the 
acceleration is reversed however maintaining the center 
of curvature. This is illustrated in Fig. [2t and[2]I. It is 
noteworthy to mention that increasing the absolute value 
of the order m results in increasing the distance between 
localized accelerating fields and center of coordinates. 
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In fact, a focused beam, associated with an incomplete 
Bessel wave field of order m = 0, can be transformed 
into an accelerating beam by a circular blazed grating. 
Furthermore, a meta-surface periodically modulating the 
amplitude of a focused wave, instead of its phase distribu¬ 
tion, will introduce a symmetrization of the accelerated 
field with respect to the y axis. The latter represents the 
basis of our study. 



FIG. 2. (a) Intensity of the magnetic field H corresponding 

to a Bessel wave field of order m — 10 propagating at A = 
632.8 nm. Accelerating wave field with (b) semi-aperture £2 = 
7t/2 and m — +10, and (c) Q — 7r/4 and m — —20. (d) 

Schematics of the circular trajectory of an incomplete Bessel 
field. 



FIG. 3. Superposition of two Bessel-driven accelerating 
beams of the same semi-aperture Q = 7r/2 and order m— 10 
and m = —10, respectively. Again A = 632.8 nm. An in- 
tereference pattern of the field intensity is evident along the 
y axis. 


III. SYMMETRIZATION OF THE INCOMPLETE 
BESSEL BEAM 

The coherent superposition of accelerating beams with 
Bessel signatures may be used to create new beam 
shapes. For instance, by combining two incomplete 
Bessel fields of the same order m but opposite sign, we 
obtain a mirror-symmetric distribution of light, as il¬ 
lustrated in Fig. [3l In this case, the phase-only term 
exp (imO) shown in Eq. ([2]) will be substituted by the 
real-valued periodic modulation 2 cos (mO). Note that an 
intereference pattern may arise in the vicinities of the y 
axis for high semi-apertures ft. 

The term exp(im$) in Eq. (|2]) behaves like a phase 
term produced by a blazed grating obtaining a maximum 
diffraction efficiency at a distance \m\/k from focal point. 


FIG. 4. Binarization of the modulation function cos (m0) and 
its influence in the Debye diffraction integral: (a) schematic 
of the binary function and (b) resulting field intensity of the 
diffracted field by using such binary function for m — 10 and 
A — 632.8 nm. 

Here we consider the generation of mirror-symmetric 
accelerating beams by using diffracting meta-surfaces 
that shape a given impinging focused wave near its focus. 
For practical reasons, the modulation of the converging 
wave will be carried out more efficiently in phase rather 
than in amplitude. Therefore we conveniently binarize 
the cosine function, as shown in Fig. 0^i. This transfor¬ 
mation will not significantly modify the main contribu¬ 
tion to the diffracted field, that is carried by the first 
diffraction order. However noisy sidelobes corresponding 
to high diffraction orders may appear. Fig. [4)3 illustrates 
the resulting field intensity \H(r)\ 2 whose phase distri¬ 
bution is governed by the binary function of periodicity 
2tt/ m. Note that the resultant intensity slightly deviates 
from that shown in Fig. [3] 


(a) (b) (c) 



FIG. 5. Intensity of the magnetic field numerically eval¬ 
uated when using a surface current (SC). In (a) we use a 
super-Gaussian apodization and in (b) and (c) we additionally 
modulate the phase distribution of the current by a piecewise 
binary function of 7r dephase. The period of the modulation 
is 2tt /m, where m — 16. 

As a proof of concept, in Fig. [5] we show some simple 
numerical experiments performed by means of a finite- 
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element-analysis commercial software^. A uniform sur¬ 
face current (SC) that is set at a given distance from 
the origin of coordinates will excite a wave field that fo¬ 
cuses at F, as shown in Fig. 0t, where the applied wave¬ 
length is 632.8 nm. By introducing a binary modula¬ 
tion of the surface current of angular period 27r/|ra|, the 
emitted field propagates near the focus exhibiting a spa¬ 
tial acceleration and a mirror symmetry with respect to 
the optical axis. The surface current distribution is ad¬ 
ditionally apodized by the super-Gaussian distribution 
A(6); we also show the response for a parameter m = 16 
at different semi-apertures: ft = 3i r/8 in Fig. 0) and 
Q = 7r/2 in (c). A periodic modulation of dephase i r 
produced in adjacent domains will produce a light cap¬ 
sule provided that the semi-aperture Q remains close to 
7r/2, that is in agreement with our analysis based on the 
Debye diffraction integral. Nevertheless, even in the op¬ 
timal case of Q = 7r/2, the hollow beam is not fully closed 
due to the proximity of the surface current to the caus¬ 
tic curve. Otherwise, increasing the radius of the surface 
current, the field distribution approaches that computed 
by Eq. (pQ). 

IV. META-SURFACE DESIGN 









^d 





(a) (b) 

FIG. 6. (a) Geometrical interpretation of the beam shap¬ 

ing using optical rays. The impinging rays split and sym¬ 
metrically bent their trajectories with the same angle. The 
resulting caustic curve is formed by two concentric circular 
segments of the same radius, (b) Design of two periodic me¬ 
dia that modulate a transmitted wave inducing a dephase of 
7r radians. For the sake of clarity we show the basic nanos- 
tructured arrangement in a planar geometry. 

As shown above, combining two incomplete Bessel 
beams with positive and negative acceleration parame¬ 
ter m may be achieved by molding the cylindrical wave- 
front of a converging wave introducing a phase-only bi- 
narization. Such beam shaping will be carried out by a 
properly-designed plasmonic meta-surface. A geometri¬ 
cal interpretation of the beam shaping using optical rays 
is illustrated in Fig. 0t. In the simple yet efficient model 
followed previously, the meta-surface will be segmented 
into regions wherein individually we induce a nearly-flat 


phase response in transmission, in a way that addition¬ 
ally neighbouring domains exhibit a dephase of it. 

In order to generate our light capsule we propose a 
locally-distributed periodic metal-dielectric media. The 
manipulation of light, such as beam acceleration or focus¬ 
ing in free space, can be achieved by controlling surface 
plasmon polaritons (SPPs) that are excited at the en¬ 
trance surface of metallic slits and corrugations with a 
subwavelength size2£r— . For instance, using metallic slit 
arrays we may control phase retardation between adja¬ 
cent slits22. However, when these nano-slits are tightly 
arranged, plasmonic modes couple leading to a (certainly 
limited) degree of homogeneization of the wave fields^. 
As a result, it can be used to induce a controlled phase 
delay. 

The dephase of wave passing through the grating is 
evaluated by p = Re(/3)d, where d is the propagation 
distance in the grating and /3 is the complex propagation 
constant computed from dispersion relation of periodic 
metal-dielectric multilayers medium for TM-polarized 
waves which is achieved by the Floquet-Bloch theorem. 
Considering a multilayered structure of period A com¬ 
posed of a metal and a dielectric of permittivity e m and 
ed, respectively, and using the transfer matrix formula¬ 
tion we finally obtain^i 

cos(iCY) = cos (ip m ) cos (<p d ) - V sin(<p m ) sin(y> d ), (4) 
where 

V = 77 -+- • ( 5 ) 

4 \^d^m tm^d J 

In Eq. (|2]), if q = KqW q , q = {m,d} refers to either the 
metal or the dielectric, 

Kq = \[tW J — /? 2 , (6) 

is the propagation constant in either the metal or the di¬ 
electric along the direction that is perpendicular to the 
metal-dielectric interface, w q denotes the width of each 
material layer, and K is the Bloch wave number. Since, in 
the simulations is considered that the incident wave im¬ 
pinges normally to the metallic grating, Eq. ® is solved 
for K = 0 (propagation all along the slits). Specifically 
the metal was gold with e m = —10.77 + 0.79^2 at the 
wavelength of interest A = 632.8 nm; also the dielectric 
under consideration was air. 

According to Eq. ®, the dephase ip gained by a Bloch 
mode is highly dependent on the width of the layers. 
Therefore, the phase retardation which is required for 
beam shaping can be obtained by controlling variation of 
the width of metal slabs and the periodicity of the array. 
Finally, an alternating sequence of two dissimilar multi¬ 
layered metamaterials, displaced transversally to the di¬ 
rection of propagation of the incident field to produce a 
beam shaping, will give rise to the required binary phase 
distribution. This is illustrated in Fig. 03. In our design, 
the width of metal is kept fixed, = 100 nm, however 
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the grating periods are Ai = 120 nm and A 2 = 142 nm. 
The separation between the entrance plane and the exit 
plane, d = 600 nm, that is the width of the meta-surface, 
remains below the wavelength. As a consequence, the 
dephase Re (A /3)d induced by adjacent photonic crystals 
will yield n radians. It should be noted that there is 
certain flexibility in order to attain the required phase 
distribution^. 

Next we will consider a subwavelength meta-surface in 
such a way that the metallic nanoelements are arranged 
over a circular sector. The effective length of the com¬ 
pound meta-surface, that is the part that will be illumi¬ 
nated by the cylindrical wave field, as measured on the 
exit surface, is equal to L = 2HR, where R is the inner 
radius of the curvilinear nanostructure. In our numerical 
simulations R = 4 fim and Q = 7r/2, and consequently the 
meta-surface has to include a minimum of m/2 bi-grating 
groups (periods). A converging super-Gaussian beam im¬ 
pinges on the circular surface of higher radius, subse¬ 
quently propagating inside metal-dielectric meta-surface 
for a distance of d = 600 nm, where the metal-dielectric 
surfaces are concentric to the focal point of the incident 
wave field. We should point out that due to the curvi¬ 
linear distribution of the meta-surface, the layer widths 
will linearly increase upon distance to the focal point F. 
However, our approach is still valid since the propagation 
constant /3 will be practically kept unaltered for a given 
metal filling fraction^. 






O 


m=28 

m=l 6 

m=10 


FIG. 8. Intensity distribution of magnetic field for different 
acceleration parameter m, but maintaining Q s= 7r/2. 


plete Bessel wave field is derived from an ideal binary 
phase modulation. To conclude we investigate the effect 
of the acceleration parameters m. Since the radius of 
Bessel beams is dependent on m, we can generate a light 
capsule with higher radius by simply increasing the value 
of m. Fig. [8] shows some numerical simulations for dif¬ 
ferent values of m, illustrating the radial tuning of the 
optical bottle. 


V. DISCUSSION AND CONCLUSION 
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Q=7r/4 

Q=3tt/8 

Cl=n/2 


FIG. 7. Intensity distribution of magnetic field for a mirror- 
symmetric incomplete Bessel beam depending on different an¬ 
gular aperture. With H — 7r/2, a nearly-closed capsule light 
can be created. 

Finally we will show the performance of the designed 
meta-surface to produce a mirror-symmetric incomplete 
Bessel beam. In Fig. [7| we show the intensity of the 
magnetic field resulting from our finite-element analysis 
for different angular semi-aperture Q. A circular (rather 
than flat) meta-surface composed of a metallic bi-grating 
is placed behind a surface current, the former inducing 
a phase modulation and the latter generating a super- 
Gaussian aberration-free focused field. We demonstrate 
that a nearly-closed capsule light can be created by us¬ 
ing a high numerical aperture focused field with H = 7r/2. 
However, the optical bottle increasingly opens for lower 
numerical apertures. 

Analyzing the results shown in Fig. [7J one can realize 
that the simulations are in good agreement with theory, 
as shown in Fig. [5] where the mirror-symmetry incom- 


In conclusion, we proposed a simple method to shape 
a high-aperture focused field, in the near field, into a 
nonparaxial accelerating field exhibiting mirror symme¬ 
try and that can be used as an optical capsule. Through 
both wave and ray optics we showed how axi-symmetric 
circular caustics can result from a binary phase modu¬ 
lation. The resultant nonparaxial hollow fields show a 
Bessel signature and have beam-widths near the diffrac¬ 
tion limit. While these fields are two-dimensional, spher¬ 
ical light shells can be easily produced by using meta¬ 
surfaces modulating the wavefront in the azimuthal co¬ 
ordinate. Also, the paths followed by the maxima will be 
circular due to the use of solutions separable in spherical 
coordinates. Our procedure can also be used for pro¬ 
ducing ID or 2D optical bottles with arbitrary convex 
caustics, which is governed by the shape of the diffacting 
meta-surface, thus providing a versatile tool for all kinds 
of accelerating waves. 

The engineered light capsules are easily sizable and ex¬ 
ert forces that can be exploited to realize micrometer-size 
optical bottles that could attract microparticles or cells. 
The beam shaping proposed here is expected to offer cer¬ 
tain advantages over the direct generation of accelerating 
beams using phase masks only, such as the ability to com¬ 
bine mechanical and optical traps produced in the vicini¬ 
ties of the focal region. Potential applications also in¬ 
clude light-induced curved plasma channels, self-bending 
electron beams, and accelerating plasmons. 
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